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Abstract 

By means of the nuclear parton distributions which can be used to 
provide a good explanation for the EMC effect in the whole x range, we 
investigate the energy loss effect in nuclear Drell-Yan process. When 
the cross section of lepton pair production is considered varying with 
the center-of-mass energy of the nucleon-nucleon collision, we find that 
the nuclear Drell-Yan(DY) ratio is suppressed due to the energy loss, 
which balances the overestimate of the DY ratio only in consideration 
of the effect of nuclear parton distributions. 

1 Introduction 

In 1983, the European Muon Collaboration(EMC) surprisingly found that 
the nucleon structure function, as measured in deep-inelastic lepton-nucleus 
scattering (DIS), varies with the target nucleus. This phenomenon has been 
known as the EMC effect. Since the discovery of the EMC effect, various 
models have been proposed to investigate the nuclear effect P, ||, Q . Among 
them, the x— rescaling P] or Q^-rescaling mechanism 0] is commonly ac- 
cepted to adquately explain the EMC effect. In addition, continuum dimuon 
production in high-energy hadron collisions, known as the Drell-Yan(DY) 
process f^, provides an independent measure of the modification of the quark 
structure of nuclei. Recently, many investigations on the EMC and nuclear 
Drell-Yan effects are still going on with great progress P, 0. Several years 
ago, the E772 Collaboration at Fermilab published the data of high-mass 
dilepton production measured in the nuclear Drell-Yan (DY) process. These 
data aroused special attention in clarifying the various explanations of the 
nuclear effect on the parton distributions. Their results show that the ratio 
of DY dimuon yield per nucleon on a nuclear target to that on a free nucleon 
(from now on, shorten as the nuclear DY ratio) is slightly less than unity 
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if the momentum fraction x carried by a target quark is less than 0.1. The 
ratios over the range 0.10 < x < 0.30, however, do not reveal distinct nuclear 
dependence. Bickerstaff et al. found that although most of the theoreti- 
cal models provide good explanations for the EMC effect, they do not give 
a good description of the nuclear Drell-Yan ratio. Most of the theoretical 
models of the EMC effect overestimate the nuclear DY ratio. 

Several years ago, we put forward an extended x— rescaling model |l^ in 
which different x— rescaling parameters for the valence quarks and sea quarks 
(gluons) in the nucleon structure function are employed in consideration of 
the nuclear momentum conservation. With a simultaneous consideration of 
the nuclear shadowing and nuclear momentum conservation, the experimen- 
tal data of the EMC effect can be well explained in the whole x region. 
However, similar to the prediction for nuclear DY ratio in pion-excess and 
quark-cluster models by using the obtained nuclear parton distributions 
in the extended x-rescaling model, the nuclear DY ratio is also overestimated. 
The difference of the nuclear effects between the nuclear DY and DIS pro- 
cesses is not clear yet. In this paper, we suggest an additional nuclear effect 
due to the energy loss in DY process. We find that the nuclear DY ratio 
is suppressed significantly as a consequence of continuous energy loss of the 
projectile nucleon to the target nucleon in their successive binary nucleon- 
nucleon collisions. This suppression balances the overestimate of the DY 
ratio only in consideration of the nuclear effect on the parton distributions. 
Therefore, a combination of these two types of nuclear effects can give a good 
explanation of the experimental data of the nuclear DY ratio. 

2 Nuclear Parton Distributions in the Ex- 
tended rescaling Model 

To provide the nuclear parton distributions which can be used to explain the 
experimental data of the EMC effect in the whole x region, we work in our fa- 
miliar extended x— rescaling model. Let Ia{n){x, Q^), I=V, S, G be the prob- 
abihty distributions of valence quarks(V), sea quarks(S), and gluons(G) in the 
nucleus A(or nucleon N), respectively. Then K^(^j^-^{x,Q^) = xIa{n){x,Q'^), 
I=V, S, G, are the momentum distributions of valence quarks(V), sea quarks(S), 
and gluons(G) in the nucleus A (or nucleon N), respectively. In Ref. [0, we 
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pointed out that, the nuclear binding effect together with the x— rescahng 
mechanism further introduced does not affect the valence quark number con- 
servation. However, the nuclear momentum is no longer conserved in the 
X— rescaling model. In order to keep nuclear momentum conservation, we ex- 
tended the x— rescaling model and employ different x— rescaling parameters 
for the momentum distributions of valence quarks and sea quarks (gluons) 
in the nucleon structure function, i.e., 

Kl^'\x,Q') = Kl^'\6y^s)X,Q'), (1) 

Because the momentum distributions of sea quarks and gluons have simi- 
lar forms, we take the same x— rescaling parameter for the momentum dis- 
tributions of sea quarks and gluons in the nucleus. The numerical result 
shows that, by properly choosing these parameters (one of them is deter- 
mined according to the nuclear momentum conservation condition), one can 
well explain the experimental data of the EMC effect. Because of its simple 
form and also for giving a good explanation of the EMC effect, the extended 



X— rescaling model has been adopted by EMC to fit their experimental 
data. Naturally, one hope that the nuclear DY ratio can also be well pre- 
dicted by using the obtained nuclear parton distributions. Unfortunately, 
the calculation results show that the nuclear DY ratio is overestimated only 
if the nuclear effect on the parton distributions is considered (see the dotted- 
line in Fig. 1). This indicates that some other mechanisms of the nuclear 
effect should be further taken into account. 



3 Energy Loss in Nuclear DY Process 

The Drell-Yan(DY) model gives a good description of the continum of 
massive dimuon pair production in the collision of proton with the nucleus 
A: 

p + A^ fi+lj,- + X. (2) 

This process is described as an electromagnetic annihilation of a quark (an- 
tiquark) in the proton p and an antiquark (quark) in the nucleon embedded 
in the nucleus A into a dimuon pair. The parton-model cross section for the 
DY process is given by 
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where i^'-factor, with ~ 2, is due to next-to leading order QCD calculations 
|l^, and a is the fine-structure constant, Cj is the fractional charge of the 



quark of fiavor i, and qi^^\x) and qi^^\x) are, respectively, the quark and 
anti-quark distributions in the proton (nucleon embedded in the nucleus A). 
The Feynman scaling variable is defined as 

xf = (4) 
y/S 

where ^/s is the nucleon-nucleon center-of-mass system(cms) energy and pi is 
the longitudinal momentum of the virtual photon of mass M. The quantities 
Xi^2 are related to Xi? and by 




xi^2 = ^{^xI + ^-±xf), (5) 

Xi — X2 = Xp, (6) 

and 

X1X2 = — . (7) 

The above relations are easily extended to account for the evolution of quark 
structure functions with Q^. 

Now let us turn to the case in which the effect of energy loss in the initial 
states is taken into account. For a nucleon- nucleus collision, the probability 
of having n collisions at an impact parameter b can be expressed as 



Pib,n) = f [T(g)a.„]"[l-T(6)a..]-^-", (8) 
nl[A — n)l 

where (Tj„ (~ 30mb |]IB[) is the non-diffractive cross section for inelastic 
nucleon-nucleon collision, and T[b) is the thickness function of the impact 
parameter b. The basic thickness function T{b) can be well approximated by 
a Gaussian function with a standard deviation Pp. If the collided nuclei are 
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small (A < 32), their density function p can also be taken to be a Gaussian 
function of the spatial coordinates. Consequently, the thickness function can 



be conveniently written as |TB 



T{h) = expi-P/2Pl)/2nPl (9) 

In terms of the standard root-mean-squared-radius parameter Tq for the nu- 
cleus A, the standard deviation Pa is given by 

Pa = r'oA'/yVs, (10) 

here Tq is found to be 1.05 fm in Ref. |jl3[, and therefore 

Pa = 0.606A^/l (11) 

For the nucleus with larger mass number A{A > 32), the thickness function 
can be approximated by using a sharp-cutoff density distribution of the form 



T{b) = ^^R\-b'e{RA-\b\), (12) 

where Ra = tqA^^^ is the radius of colliding nucleus with ro = 1.2 fm. 

In Eq. (|^), the first factor on the right-hand side represents the number 
of combinations for finding n collisions out of A possible nucleon-nucleon 
encounters, the second factor gives the probability of exactly n collisions and 
the third factor gives the probability of having exactly A — n misses. The 
total probability for the occurrence of an inelastic event in the collision of 
proton with the nucleus A at an impact parameter b is the sum of Eq. (|^) 
from n = 1 to n = ^: 



= Y.P{n,b) 

db n=l 

A 



= 1-[1-T(6)ai„]^. (13) 

Therefore, from Eq. (|13D, the total inelastic cross section cr^n^ for the colli- 
sions of protons with the nucleus A is 

c^fn^ = j db{l-[l-T{b)a,,,f}. (14) 
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In an inelastic nucleon-nucleus collision without impact parameter selection, 
the number of nucleon-nucleon collisions n (for n = 1 to A)has a probability 
distribution P{n). This is obtained by integrating P{n,b) over all impact 
parameters: 

Pin) = JdbPjnrb) (^5^ 
E JdbP{n,b) 

n=l 

where the denominator is to ensure that P{n) is properly normalized as 

E Pin) = 1. (16) 

n=l 

From Eq. (|T^), the denominator of the right hand side of Eq. (|T3p can be 
replaced by 

E / dbPin, b)= [db{l-[l- T(6)a,„]^}. (17) 

n=l 

To describe the energy loss in the collision of a proton with the nucleus A, 
we start with remarks on the relative role of "soft" and "hard" interactions 
in nuclear collisions at very high energy: The incident proton interacts with 
spectator nucleon and makes soft (nonperturbative) minimum bias collisions 
before making the high dimuon pair; During the "soft" collisions, the 
projectile proton imparts energy to the struck nucleon and therefore must 
loose energy; Thus energy loss must affect the cross sections of producing 
dimuon pair. After the projectile proton has had additional n collisions with 
nucleons embedded in the nucleus, the cms energy of the colliding nucleons 
with "hard" DY collisions can be expressed as 

where generally taken as 0.2-OAGeV , is the cms energy loss per collision 
in the initial state. Therefore, the cross section for the DY process can be 
re-written as 

d'a 47ra^ ^ , [q!ix[)qfix',) + qKx[)qfix',)] 



dM^dxp Vs' 9s' i * Jx'Z + AMys' 
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where the rescaled quantities are defined as 



x'p = —tL = rsXF, (20) 



and 



with the cms energy ratio: 



x'l 2 = r^xi,2, (21) 



(22) 



The average cross section for the dimuon production in nuclear DY process 
can be expressed as 

d'^a ^ (fa 

) = i:Pin)J^_. (23) 



To make a comparison between the theoretical prediction for the nuclear DY 
ratio and the experimental data with respect to the variables xi and X2, Eq. 
(p!9D and Eq. (ESf) can be re-expressed as 



"^"^ Ktu^ll^M{rsXi)qf{rsX2) + (fi{rsXi)qt{r,X2)l (24) 



dxidx2 9M2 

and 



d^a ^ d^a 

:^) = T.Pin)-^, (25) 
dxidx2 ~i dxidx2 



respectively. It is noteworthy that in ( p2D is always greater than one if 
there exists an energy loss in the collisions of protons with the nucleus A. 



4 Numerical Results 



To compare the theoretical prediction for nuclear Drell-Yan process with the 
experimental data of E772 collaboration f^, we introduce the nuclear DY 
ratio as 

J dX\ 'fln^-^fl^^ 
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where ^ is the differential cross section for the dimuon pair production 



in the proton-deuteron colhsion. The integral range for xi in Eq.(26) is 
determined according to the kinematic region of the experiment in Ref. [^], 
i.e. X1—X2 > 0, and 0.025 < X2 < 0.30. By using the free parton distributions 



of GRV and taking the x— rescaling parameters Sy = 1.026 and Ss = 
0.945 for valence quarks and sea quarks, respectively, we calculate the nuclear 
Drell-Yan ratios for ^^Fe. In addition, we assume that the gluon and the 
sea quark structure functions have similar shadowing effects at small X2, and 
according to Ref. [l^ , we introduce the shadowing factor i?^^ both for gluons 
and sea quarks as 

A / N _ / 1 + alnAln{x2/ 0.08), X2 < 0.08, 

^sh[X2) - j 1 + blnAln{x2/0.08)ln{x2/0.2A), 0.08 < X2 < 0.3 

where the parameters a, b are taken as 0.025 and -0.02, respectively. The 
calculation results with = 0.0, 0.2, O.AGeV shown in Fig. 1 indicate that 
the energy loss effect is essential in the explanation of the nuclear DY ratio. 



5 Discussion and Summary 

To sum up, by means of the nuclear parton distributions in the extended 
x— rescaling model on which the experimental data of the EMC effect can 
be well explained, we investigate the nuclear DY process focusing on the 
continuous energy loss of the projectile nucleon to target nucleons in their 
successive nucleon-nucleon collisions. We find that the nuclear DY ratio will 
be overestimated if the nuclear effect on the parton distributions is the only 
factor considered. The calculation results show that the nuclear DY ratio is 
sensitive to the change of sea quark distributions. On one hand, to adequately 
explain the antishadowing effect and the EMC effect, there should exist an 
enhancement of the sea quarks in the nuclear quark distributions in the range 
of 0.1 < X < 0.3, resulting in the overestimation of the nuclear DY ratio; On 
the other hand, the energy loss causes a suppression of the nuclear DY ratio, 
which balances the overestimate of the nuclear DY ratio due to the nuclear 
effect on the parton distributions. So, in the nuclear DY process, there is no 
distinct nuclear effect as observed in DIS process due to the the combination 
of the nuclear effect on the parton distributions and the energy loss effect. 
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Similarly, the J/'i/; suppression is also partially due to the energy loss in the 
initial states. 
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Figure Caption 

Fig. 1. The nuclear Drell-Yan ratios T^^^^{x2) predicted with the energy 
loss ( solid curve for ^^=0.2 GeV, dashed curve for ^^=0.4 GeV) and 

without the energy loss(dotted curve for ^^=0.0 GeV). The experimental 
data are taken from the E772 Collaboration [Q. 
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